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ABSTRACT 

Episodic activity of quasars is driving growth of supermassive black holes (SMBHs) via accretion of baryon 
gas. In this Letter, we develop a simple method to analyse the duty cycle of quasars up to redshift z ~ 6 universe 
from luminosity functions (LFs). We find that the duty cycle below redshift z ~ 2 follows the cosmic history of 
star formation rate (SFR) density. Beyond z ~ 2, the evolutionary trends of the duty cycle are just opposite to that 
of the cosmic SFR density history, implying the role of feedback from black hole activity. With the duty cycle, we 
get the net lifetime of quasars (z < 5) about ~ 10 9 yrs. Based on the local SMBHs, the mean mass of SMBHs is 
obtained at any redshifts and their seeds are of 10 5 M Q at the reionization epoch (z re ) of the universe through the 
conservation of the black hole number density in comoving frame. We find that primordial black holes (~ 1O 3 M ) 
are able to grow up to the seeds via a moderate super-Eddington accretion of ~ 30 times of the critical rate from 
z = 24 to z re . Highly super-Eddington accretion onto the primordials is not necessary. 
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Subject headings: black hole physics - 
general 

1. INTRODUCTION 



galaxies: active — galaxies: evolution — galaxies: nuclei — quasar: 



Accretion of gas onto SMBHs is powering the huge energy 
output from quasars, but SMBH formation and quasar's life- 
time remain open in this well-established paradigm. The el- 
egant idea, comparing the total accreted mass density during 
the active phases with the local mass density of SMBHs in nor- 
mal galaxies (Soltan 1982, hereafter Soltan argument), has been 
examined in detail from quasar LF (Chokshi & Turner 1992; 
Yu & Tremaine 2002) and X-ray background (Marconi et al. 
2004). This generally convinces us that accretion during their 
episodic activities is the main source of mass growth. However, 
the SMBH growth at different redshifts remains open. 

In the present paper, we make an attempt to develop an effi- 
cient way to estimate the duty cycle for high redshifts, offering 
a new clue to understand SMBH growth at different redshifts. 
We use the cosmological parameters Hq = 75 km s" 1 Mpc -1 , 
f^M = 0.3 and J7a = 0.7 throughout this paper. 

2. DUTY CYCLE OF QUASARS 

The SMBH mass function at redshift z in active and dormant 
galaxies is J\f(M m , z), the duty cycle is defined as a fraction of 
active black holes to the total, 



S(M.,z) = 



(1) 



M{M.,zY 

where M, is black hole mass and ^>{M,,z) is the mass func- 
tion of black holes in quasars at a redshift z. The averaged duty 
cycle is 8\(z) = N qso (z)/N a }i(z) in term of the number density, 
where N qso (z) = J M , *(M.,z)dM„ jV a „(z) = J M . A/"(M.,z)dM„ 
and M* is the lower mass limit of SMBHs in the sample. 
In the meanwhile, the mass-weighted duty cycle is given by 
Si(z) = P.{z)/ P f{z\ where p.(z) = J M , *(M.,z)M.dM. and 

pf( z ) = j M ,J\f(M.,z)M,AM.. It has been shown that 6i(z) = 
#2(z) = S(z) in Wang et al. (2006c), namely 

P»(z) 



S(z)-- 



Pt(z) 



(2) 



Eq. (2) involves that the mean masses of active and inactive 
black holes are equal, and then converts the number density to 
mass density ratio (Wang et al. 2006c). It should be noted 
that <5(z) represents the duty cycle of the major population of 
SMBHs at a given redshift. 

We get the mass density of black holes from quasar LF if we 
assume a constant in of quasars as in literatures (e.g. Marconi 
et al. 2004, but see evidence for this in Kollmeier et al. 2006). 
The dimensionless accretion rate is defined as m = M/M a - lU 
where M is accretion rates of black holes, M cr ; t = L^/c 2 the 
critical rate, Lead = M \c 2 / 'fs a i p the Eddington luminosity, ?s a i P = 
arc I '4TrGm p = 0.45 Gyr the Salpeter time, <tj the Thomson scat- 
tering section, c the light speed, G the gravity constant and m p 
the proton mass. With the help of LboI = i]Mc 2 = r\rhM % c 2 /fs a i P , 
where 77 is the radiative efficiency and quasar LF, we have the 
black hole mass density in active galaxies at redshift z 



P.(z)-- 



1 



rj{m)c- 



■U(z)t Sa i v , 



(3) 



where (m) is the mean dimensionless accretion rate, and the 
luminosity density is given by 



U(Z)= / LBol$(iBol,z)dL B ol, 



(4) 



where $(Lbo1,z) is the bolometric LF, LboI the bolometric lumi- 
nosity, and Lg ol its corresponding limit. If the mass density of 
seed black holes of quasars is pi at their birth epoch (z max ), the 
mass density of all (active and inactive) black holes is given by 



Pf(z)-- 



P S .+ 



l-yU(z) 
77 c 2 



pI + 



\-r,U{z) 
77 c z 



(5) 



where U(z) = J z <mi " U(dt /dz)dz. The dependence of 77 on z can 
be neglected in eq. (5) as shown by Wang et al. (2006b) from 
Sloan Digital Sky Survey (SDSS) data. We almost know noth- 
ing about pi at z max except for some limited information on the 
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reionization of the universe (Madau et al. 2004). Inserting (3) 
and (5) into (2), we have 

, _ t/(z)?Salp _ U (z)fSalp 

(Z) ~ (m) [Us + ( 1 - rj)U(z)] ~ (m) ( 1 - V ) U(z) ' 
where Us = ijplc 2 is the energy density of the seed black holes 
and the approximation is valid for p\j ' p^ l {z) <§C 1. 

Information on p\ can be estimated from the reionization of 
the universe. The WMAP (Wilkinson Microwave Anisotropy 
Probe) detected a large optical depth to Thomson scattering, 
r e = 0.17 ±0.04, and suggests that the reionization happened at 
much higher redshift z re = 17 ± 3 (Spergel et al. 2003). Madau 
et al. (2004) suggest that the reionization may be powered 
by mini-quasars, which have a mass density of accreting black 
holes at least 2 x 10 3 M Q /Mpc 3 , which roughly agrees with that 
extrapolated by the luminosity function used below. We take 
pl=2x 10 3 M Q /Mpc 3 and 77 = 0.1 in this paper. 

We would like to point out followings: 1) the approxima- 
tion is accurate enough within a certain redshift z e when pi can 
be neglected and baryon accretion dominates. It breaks when 
S > 1 ; 2) the advantage of eq. (6) is that S is not sensitive to 77; 
3) 5 is insensitive, if the specific LF is applied, to the bolomet- 
ric correction factor since it will be canceled on both sides of 
the numerator and the denominator; 4) (m) is an observable in 
principle and seems to be a constant at least between z = 0.3 ~ 4 
(Kollmeier et al. 2006, but see Netzer et al. 2007 for a small 
high redshift sample). Though the cosmological evolution of m 
is poorly understood (Netzer & Trakhtenbrot 2007), the influ- 
ence of (m) is clear in eq. (6). These make eq. (6) robust to 
calculate <5 for high redshift SMBHs and accurate enough for 
low redshift ones. 

In this Letter, we use the bolometric LF given by Hopkins 
et al. (2007). It is combined through bolometric luminosity 
correction from a large set of LFs in optical, soft and hard X- 
rays, and near- and middle-IR bands and also covers the fraction 
of obscured quasars (see also Maiolino et al. 2007; Miiller & 
Hasinger 2007 for the fraction of type II quasars). We use the 
LDDE LF given by eq. (1 1-16), of which parameters are listed 
in Table 4 in Hopkins et al. (2007) and extrapolate luminosity 
functions beyond z = 6.0. We assume z max = z re and find that the 
final results are not sensitive to z max . 

Fig la shows t/(z). There is a break at z ~ 2 (also small 
effects on the duty cycle), which is caused by the LF. U dra- 
matically drops toward low redshifts and gradually decreases 
toward high redshifts. The function U(z) is not plotted here, but 
its trends are equivalent to that of the cumulative mass density 
of black holes, whose behaviors can be seen from (M,(z)) ac- 
cording to eq. (9) (shown in Fig. 2b). Quasars with redshifts 
of 0.3 < z < 4 have a mean Eddington ratio of LBoi/^Edd = 0.25 
with a scatter of 0.3 dex, namely, (m) = 2.5t]q\, where 770.1 = 
77/O.I, using the empirical reverberation relation (Kollmeier et 
al. 2006). Fig lb shows the duty cycle of quasars ranging from 
z = to 6 for (m) = 1,2.5,5, 10, respectively. 

From Fig. lb, we find that there is a gradual decrease of 5 
from high redshifts to z ~ 2 and then dramatically evolves go- 
ing down to 10~ 3 ~ 10 -4 in the local universe. The duty cycle 
below z < 2 is similar to the results from Wang et al. (2006c). 
Comparing with the cosmic history of SFR density, we find 5 
has very similar evolutionary trend. Exhaustion of gas leads to 
lower psfr, so does the fueling gas to SMBHs. This strong evo- 
lution is thus regarded as the lack of fueling gas, consequently, 
most SMBHs are then starving. Afterglows of quasars are ex- 
pected to appear (Wang et al. 2005). 




redshift (z) 

FIG. 1 . — (a) The luminosity density versus redshift based on bolometric 
luminosity function (Hopkins et al. 2007). {b) The duty cycle versus redshifts. 
We calculate cases of (m) = 10,5,2.5, 1, corresponding to that quasars are ra- 
diating at LEdd, 0.5LeiM> 0.25Lrj,id and 0. lZ-Edd, respectively. The error bars are 
taken from the averaged value of LF A<5/<$ = A<I>/<I> fa 0.2. The cosmic his- 
tory of SFR density is inserted as open circles taken from Reddy et al. (2007), 
but scaled by a factor of 10 3 Mq yr~' Mpc -3 . 

Beyond redshift z ~ 2, the trends of <5 and psfr are just op- 
posite as shown in Fig lb. The present S at high redshifts agree 
very well with the duty cycle estimated from clustering of high 
redshift quasars (z > 2.9) by Shen et al. (2006). The duty cy- 
cle gradually decreases from high to low redshifts while psfr 
slowly increases. During this epoch there is enough gas for both 
star formation and accretion onto SMBHs. We note that the 
strong feedback of black hole activities (Schawinski et al. 2006; 
Wang et al. 2007) is driving to blow away fueling gas so as to 
switch off quasars (Di Matteo et al. 2005). It has been found by 
Peng et al. (2006) and McLure et al. (2006) that quasars at z - 2 
have stellar mass less than expected from the local Magorrian 
relation. This also is indicated by the broken relation between 
star formation and AGN activity at high luminosities (Maiolino 
et al. 2007), suggesting the SMBHs are growing faster than 
star formation at high redshifts in the presence of strong feed- 
back of black hole activity. On the other hand, one can check 
if the star formation still obeys the Kennicutt-Schmidt's law as 
done in Seyfert galaxies (Wang et al. 2007). It is expected that 
further evidence for the evolution of the feedback is to search 
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from statistics at high redshifts. Future high spatial resolution 
observations of ALMA (Atacama Large Millimeter Array) will 
directly uncover the detailed nature of feedback. 

3. ACCRETION AND GROWTH 

3.1. Lifetime of Quasars 

The net lifetime of quasars (i.e. the total over the Hubble 
time) is given by 



fQSO : 



(7) 



where t is the cosmic time. Fig. 2a shows the net lifetimes 
for quasars with different accretion rates. For the rate m = 2.5, 
we find the typical value is < 10 9 yrs for z < 5. Marconi et al. 
(2004) obtained a lifetime of a few 10 8 yrs, based on m = 10, 
which agrees with our results for the same m. If a single 
episodic lifetime is measured from the transverse proximity ef- 
fect (Goncalves et al. 2007), the cycles of episodic activities 
can be estimated. The cutoff around z ~ 15 in Fig. 2a is caused 
by setting the birth of quasars at this epoch. 

3.2. Growth of Seed Black Holes from the Primordial 

Neglecting mergers, we have the conserved number density 
in the co-moving frame as 



p?{z) = P l P l 

(M.(z)) (MS) <M°) ' 



(8) 



where p° % and (A/2) are the SMBH mass density and mean 
mass in the local universe. This is justified by the results 
from detailed numerical simulations, which show the mass con- 
tributed from major mergers is only roughly a few percent 
after z = 10 epoch (Volonteri et al. 2003). The mass den- 
sity of the local SMBHs is given by p°, = 4.2 x 10 5 M Q /Mpc 3 
(Shankar et al. 2004; Marconi et al. 2004). We convert the 
function of dispersion velocity [<fi((r)] of early type galaxies 
(Sheth et al. 2003) into the mass function of the SMBHs, 
<f>(M.) = (j)(a)da/dM„ where M./M = 10 813 (cr/200) 402 is 
used (Tremaine et al. 2002). The mean mass of local SMBHs is 
given by (A/2) = / M.<j>(M.)dM./ J 4>{M.)6M. = 8.9 x 1O 7 M . 
We thus have the mean mass of black holes at z 



(M.(Z)) : 



pT(z) 



(A/2). 



(9) 



Fig. 2b shows (M,(z)) as a function of redshifts. If future ob- 
servations could provide the mean mass of SMBH [(M° bs (z))] 
for a complete sample at a given redshift, it becomes feasible 
to justify if the non-baryon accretion onto the black holes is 
necessary by a simple comparison of (M* s (z)) and (A/.(z)). 
From eq. (8), we have the mean mass of seed black holes 



{Ml)=[^)(M .)^2.0xl0 5 M Q , 



(10) 



where we use pl limited by the WMAP. It should be noted that 
the real (M^ ) could be smaller than the given by eq. (10) if stars 
partially ionize the universe. The value of (Ml) agrees well 
with the black hole mass of the mini-quasar model in Madau et 
al. (2004). How to form such a massive seed black hole at z re 
remains open from the primordial at z ~ 24. 
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FIG. 2. — (a) The net lifetime of quasars, (b) The mean mass of supermassive 
black holes at different redshifts. 

The primordial black holes can be generally produced by the 
collapse of population III stars (Madau et al. 2001), or col- 
lapse of primordial gas clouds (Haehnelt & Rees 1993; Loeb 
et al. 1994) or gravitational core collapse of relativistic star 
clusters (Volonteri 2006) with a typical mass of 10 2 - 10 3 M Q , 
10 3 - 10 6 M Q and 10 2 - 10 4 M Q , respectively. For a growing pri- 
mordial black hole with a typical mass of (Ml) = 10 3 M Q , the 
necessary accretion rate is 



fSalp 



(l-?7) A* 



In 



(Ml) 
(Ml) 



;30, 



(11) 



in the Salpeter growth, where Af (« 0.09Gyrs) is the interval 
between z = 24 ~ 17. This is only moderate super-Eddington 
and is realistic in high redshift universe. The photon trapping 
effects make the accretion have low radiative efficiency, but 
disks of the black holes radiate at a level of L Edd (Wang et al. 
1999; Wang & Zhou 1999; Ohsuga et al. 2005). 

We note that this conclusion is inconsistent with that in 
Volonteri & Rees (2005, hereafter VR05). VR05 suggested that 
a highly super-Eddington accretion onto primordial black holes 
between z = 23 <~ 24 to form seed black holes of 10 5 M Q in term 
of the Bondi accretion (see their Fig. 1), and subsequently the 
seeds gradually grow up to be a 10 9 M Q SMBH until z = 6. 
Using the Bondi rate, we obtain m = m Q m l l/ [l-mj (f/ T o)]> 
where m® is the initial mass of a primordial black hole in 

GM Q ccr T no/c 3 = 0.05 T^ 2 ti4 and 



units of solar mass, mo 
to = c 3 /47rG 2 m p noA/© = 13.9 r g «4Gyr, n\ 
the number density of ambient medium, c s = (kT/m v )''" the 



no/10 4 cm 3 is 

1/2 



sound speed of the medium, To.s = T /8000K the temperature of 
the medium and k the Boltzmann constant. The rate m goes 
to 10 4 " 5 , even infinity for a critical time of t c = To/m®, be- 
fore the Bondi radius is larger than the typical dimension of 
the primordial clouds (VR05). The classical Bondi accretion is 
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valid provided the growth timescale of the accreting black hole, 
tbh = (dlnM./df) is much longer than the gaseous dynam- 
ical Tdy n = RBondi/cs, where ^Bondi is the Bondi radius. Under 
the environment in high redshift universe, we have T BH /7d yn = 
c 6 s /4TrG 3 M 2 m m p n = 0.1 T^Mfrq 1 , where M 5 = M./10 5 M Q . 
This directly indicates that the Bondi approximation is broken 
for a rapid growth of the primordials. A self-consistent treat- 
ment of time-dependent Bondi accretion onto a growing black 
hole is needed for such a context. 

On the other hand, the rapid growth through super-Eddington 
accretion with the Bondi rates in VR05 is suffering from 
strong feedback from outflow or Compton heating (Wang et al. 
2006a). The strength of feedback depends on the angular mo- 
mentum of the primordial gas accreted onto the black hole. A 
potentially efficient way to avoid the strong feedback is that the 
primordial gas has angular momentum large enough to extend 
the outer radius of the accretion disk so as to suppress the feed- 
back due to the Compton heating. A moderate super-Eddington 
accretion is thus feasible to form a seed black hole from the 
primordial during the period of z = 24 to 17. 



4. CONCLUSIONS AND DISCUSSIONS 

We develop a convenient way to calculate the duty cycle 
based on LF, which applies to any redshifts. We find that the 
trends of the duty cycle and the cosmic history of the SFR den- 
sity are just opposite when z > 2. This could be explained by 
AGN feedback to star formation. With the duty cycle, the net 
lifetime of quasars can be obtained from the quasar luminos- 
ity function. The mean mass of the seed black holes is up to 
2 x 10 5 Mq at z ~ 17, which is able to grow up from the pri- 
mordial (~ 10 3 Mq) at z = 24 via moderate super-Eddington ac- 
cretion. More deeper surveys are expected to improve the LFs 
for more sophisticated investigations of growth of black holes 
in high redshift universe. 
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